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摘 要 : 在 建立 雷达 高 度 计 海 况 偏差 (Sea State Bias，SSB) 非 参数 模型 时 ， 通 常会 用 到 局 部 线性 回归 (Local Linear 
Regression，LLR) 估 计 器 ， 而 传统 的 局 部 线性 回归 估计 器 涉及 高 维和 矩阵 运算 ， 当 建 模 的 数据 量 较 大 时 ， 估 计 海 况 偏 
差 需要 大 量 的 时 间 ， 从 而 使 得 非 参 数 估计 方法 很 难 用 于 高 维 海 况 偏差 模型 。 该 文 提 出 一 种 改进 的 局 部 线性 回归 
(Improved Local Linear Regression, ILLR) 估 计 器 ,可 以 避免 传统 的 LLR 估计 器 所 需 的 高 维和 矩阵 运算 ,在 不 影响 海 况 
偏差 估计 结果 的 条 件 下 ， 将 局 部 线性 回归 估计 器 获取 加 权 函 数 的 时 间 复 杂 度 由 OCN”) 降低 为 O(N) ， 从 而 大 幅 地 
降低 估计 海 况 偏差 所 需 的 时 间 ， 为 实现 高 维 非 参 数 海 况 偏差 模型 的 实时 运算 莫 定 了 基础 。 
关键 词 : 雷达 高 度 计 ; 海 况 偏差 ， 非 参数 估计 ; LLR 估计 器 
中 图 分 类 号 : TN953 文献 标识 码 : A 文章 编号 : 
DOI: 10.11999/JEIT151280 


mp 


Pee 


Improved Local Linear Regression Estimator and Its Application to the 
Estimation for Radar Altimeter Sea State Bias 


JIANG Maofei®””®” XUKe®” LIUYaong” WANG Lei®® 
VKey Laboratory of Microwave Remote Sensing, Chinese Academy of Sciences, Bewing 100190, China) 
C J National Space Science Center, Chinese Academy of Sciences, Beijing 100190, China) 
ed | University of Chinese Academy of Sciences, Bewing 100049, China) 


(Tantat Marine Environmental Monitoring Center Station, State Oceanic Administration, Yantai 264000, China) 


Abstract: The Local Linear Regression (LLR) estimator is usually used when developing a nonparametric model for 
radar altimeter Sea State Bias (SSB). However, the conventional LLR estimator contains matrices with high dimension. 


When a large number of data are used to develop the SSB model, the SSB estimation costs too much time. Therefore, the 


nonparametric estimation method can hardly be used to develop high-dimensional SSB models. This paper presents an 

_ improved LLR estimator, which can avoid high-dimensional matrix operations. The improved LLR estimator can greatly 

( reduce the time for SSB estimation without affectng the estimated accuracy. So the improved LLR estimator can laid the 
foundation for high-dimensional SSB models. 
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1 引言 


全 球 海平 面 上 升 己 越 来 越 引起 人 们 的 关注 ， 雷 达 高 度 计 的 一 个 重要 应 用 就 是 测量 平均 海面 高 度 [, 写 。 但 
是 由 于 海 况 偏 差 的 存在 ， 高 度 计 测 量 得 到 的 平均 海面 会 低 于 实际 的 平均 海面 ， 必 须 对 它 进行 校正 。 海 况 偏 
差 主要 是 由 于 海面 的 非 高 斯 分 布 特性 引起 ， 包 括 电磁 偏差 和 偏 斜 度 偏差 9。 利 用 现 有 方法 得 到 的 海 况 偏差 
的 不 确定 度 能 够 达到 2 cmr9， 已 经 成 为 Jason 系列 测 高 卫星 最 大 的 误差 源 尼 外。 

目前 ， 业 务 化 运行 的 雷达 高 度 计 的 海 况 偏差 校正 主要 采用 以 高 度 计 测 得 的 风速 和 有 效 波 高 作为 输入 的 
非 参 数 SSB 模型 [6"]。 非 参数 SSB 模型 的 构建 通常 需要 用 到 局 部 线性 回归 估计 器 ,而 传统 的 局 部 线性 回归 估 
计 器 需要 大 量 的 矩阵 运算 ， 从 而 消耗 大 量 的 运算 时 间 ， 极 大 地 限制 了 非 参数 估计 在 高 维 SSB 模型 中 的 应 用 。 

本 文 提 出 了 一 种 改进 的 局 部 线性 回归 估计 器 ， 该 估计 器 避免 了 高 维和 矩阵 运算 ， 将 从 局 部 线性 回归 佑 计 
器 获取 加 权 函 数 的 时 间 复 杂 度 由 O(N ”) 降低 为 OOV) 。 数 值 实验 结果 表明 ， 在 其 它 条 件 都 相同 的 条 件 下 ， 和 
传统 的 局 部 线性 回归 估计 器 相 比 ， 改 进 的 局 部 线性 回归 估计 器 不 会 影响 SSB 估计 的 结果 ， 但 可 以 大 幅 地 减 
低估 计 SSB 所 需 的 时 间 。 
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2 局 部 线性 回归 估计 器 及 其 在 海 况 偏差 非 参 数 估计 中 的 应 用 

GASPAR 等 人 [9] 在 1998 年 首次 将 非 参 数 估计 的 方法 用 于 SSB 模型 的 构建 ， 和 参数 模型 以 9 相 比 ， 非 参 
数 模型 的 精度 更 高 ， 从 而 被 更 广泛 地 应 用 于 高 度 计 测量 海面 高 度 的 海 况 偏差 校正 。 用 于 非 参数 SSB 模型 的 
估计 器 主要 有 Nadaraya-Watson (NW) 估 计 器 [I 和 局 部 线性 回归 估计 器 09, NW 估计 器 在 接近 边界 处 会 出 现 
较 大 偏差 ， 而 局 部 线性 回归 估计 器 则 可 以 很 好 解决 这 一 问题 ww， 所 以 局 部 线性 回归 估计 器 得 到 了 更 广泛 
的 使 用 。 
2.1 局 部 线性 回归 估计 器 

给 定 N 组 观测 值 (,x,)，x, 是 含有 p 个 变量 的 向 量 ，x, 和 y 满足 : 
y,= r(x,)+ e (1) 


中 是 回归 函数 ，e, 是 一 个 零 均 值 的 误差 项 。 对 于 回归 函数 r(x)， 在 x 的 邻 域 z 对 其 进行 Taylor 展开 有 : 
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r(z)= DT+ 3 b(z,- Xx,) (2) 


j=1 


对 于 给 定 的 N 组 观测 值 0,x,) ,可 以 通过 加 权 最 小 二 乘 的 方法 找到 一 组 5 = (5b,,b,,…,5b,)”， 使 得 式 (3) 最 小 : 


Db,- 加 -Co x)TK, Ce,- x) (9) 


式 中 ，K, (x, - x) 是 加 权 函 数 ， 对 于 局 部 线性 ] 核 函数 ， 常 用 的 核 函 数 有 高 


斯 核 函 数 和 Epanechnikov 核 函 数 。 令 y = Doyv]， 


W = diag{K,(x,- x)} 是 一 个 (W” WN) 的 加 权 和 矩阵 ， 式 (3) 可 以 改写 为 : 


0 -Xb) WO -Xb) (4) 
最 小 化 式 (4)， 可 以 得 到 : 
此 = (WI, XWy (5) 
r(x) 的 局 部 线性 下 归 估 计 器 Buntx) 为 此 ， 而 太 的 其 它 成 分 是 函数 7 的 一 阶 偏 导 数 ， 即 
Pig)= = eT WX, XWy (6) 


式 中 ,e 是 一 个 含有 p+1 个 变量 的 单位 向 量 ，e， = [1,0,.…,0]。 


对 于 给 定 的 N 组 观测 值 ,x,)， 根 据 非 参 数 估计 理论 ， 条 件 期 望 E[y |x] 表 示 在 给 定 x 的 条 件 下 y 的 均 
值 ， 其 核 估计 器 jx) 可 以 看 作 是 观测 值 y 的 加 权 平 均 : 


N 


Fo) = & yavG -zi) (7) 


生 回归 估计 器 Bur(x)， 在 式 (6) 中 ， 令 7 = eT(XJWX,) 对 W ， 就 可 以 得 到 权 函 数 : 
avy(x- x,)=T(0) (8) 
2.2 局 部 线性 回归 估计 器 在 海 况 偏差 非 参 数 估计 中 的 应 用 
通常 使 用 交叉 点 处 的 海面 高 度 不 符 值 来 建立 SSB 模型 ,卫星 绕 地 球 运行 一 周 的 轨迹 包括 两 个 弧 段 (Pass)， 
分 别 为 上 升 轨道 和 下 降 轨 道 ， 上 升 轨道 和 下 降 轨 道 的 交点 就 称 为 交叉 点 ， 而 交叉 点 海面 高 度 不 符 值 是 指 上 
升 轨道 和 下 降 轨道 在 交叉 点 处 测 得 的 海面 高 度 之 差 [13]。 
SSH 表示 未 经 SSB 校正 的 海面 高 度 测量 值 ， 假 设 SSB 是 输入 变量 x = (U,SWH) 的 函数 i ， 并 令 
y = SSH, - SSH, ， 根 据 文献 [3] 和 文献 [1o] 可 知 : 
y=j(x,)- j(x)+e (9) 
式 中 ，e 是 一 个 零 均 值 的 误差 项 ， 角 标 1 和 角 标 2 分 别 对 应 上 升 轨道 和 下 降 轨 道 。 在 给 定 x, = x 的 条 件 下 ， 


A 


对 于 局 部 线 ' 


y 的 条 件 期 望 为 : 
ED |x,=x]=j(x)- EDV (x)|x,= x] (10) 
对 于 N 组 测量 值 Q,wx,,x,)， 利 用 式 (7) 和 式 (10) 可 以 得 到 : 


一 


j 6)= ja x)+ jy- x,) (0) 


式 (11) 给 出 了 估计 j (x) 的 方法 ， 只 要 知道 了 j Ce ) ， 就 可 以 估计 出 任意 输入 x 所 对 应 的 SSB 估计 值 cc) 。 利 
用 x 谷 换 式 (11) 中 x ， 式 (11) 可 以 重 写 为 : 


将 式 (12) 转 化 为 矩阵 的 形式 可 以 得 到 ; 
-Ay ,=Ay (13) 
式 中 ， 了 和 4 都 是 (W “NN) 的 和 矩阵， 矩阵 T 是 一 个 单位 矩阵 ， 而 矩阵 4 的 元 素 为 a = a - x,)， 
jo = GD CE) y= Dy ye 
和 矩阵 I - 4 是 一 个 奇异 矩阵 ， 其 秩 为 N-1。 为 了 消除 不 确定 性 ， 可 以 先 给 7 , 中 的 第 一 个 元 素 一 个 假设 
值 : j (x, )=j。， 那 么 式 (13) 可 以 重 写 为 : 


BJ =Ay- By, (14) 
式 中 , j 包含 ji 除 j Cj) 之 外 的 N-1 个 元 素 , 即 j ”= (x,),…,j (rw)]，B, 是 算 阵 I -4 的 第 1 列 ，B| 是 
其 余 各 列 组 成 的 矩阵 ， 利 用 最 小 二 乘法 求解 7 : 
= (BBC -Bo no) (15) 
再 联合 j (x,,) ， 就 可 以 得 到 j , ， 将 它 带 入 式 (11) 就 可 以 得 到 任意 x 所 对 应 的 j () 的 估计 值 。 
在 估计 SSB 时 ， 可 制作 一 个 SSB 估计 值 关 于 风速 U、 有 效 波 高 SWH 的 2 维 查 找 表 ， 再 通过 插值 的 方 
法 得 到 任意 高 度 计 测量 点 的 SSB 估计 值 . 为 了 得 到 SSB 估计 值 的 查找 表 , 通常 是 先 利用 一 个 cycle 的 数据 计 
算 查找 表 中 每 个 网 格 点 的 输入 变量 x 所 对 应 的 SSB 估计 值 (x)， 再 对 多 个 cycle 得 到 的 结果 进行 平均 。 每 个 
cycle 的 交叉 点 的 有 效 个 数 一 般 可 达 5000~7000 个 ， 如 果 将 所 有 交叉 点 都 用 于 建 模 ， 那 么 式 (15) 中 矩阵 求 道 
需要 大 量 的 运算 ， 消 耗 大 量 计算 时 间 。 在 文献 [8] 的 模型 中 ， 从 每 个 cycle 随机 抽取 500 个 有 效 的 交叉 点 来 
建立 SSB 估计 模型 。 为 了 避免 了 式 (15) 中 算 阵 求 逆 运 算 ，GASPAR 等 人 oo 在 2002 年 利用 LSQR 算法 来 求 线 
性 方程 式 (14) 的 解 ， 并 将 每 个 cycle 的 所 有 交叉 点 数据 都 用 于 SSB 的 建 模 。LSQR 算法 是 PAIGE 等 人 [9 在 
1982 年 提出 的 一 种 解 稀 玻 线性 方程 组 的 欠 代 方法 , 它 被 用 于 求 线性 方程 组 Ax = 5b 的 解 或 者 使 [Ax bl 最 小 
的 解 。 为 了 进一步 提高 求解 线性 方程 式 (14) 的 速度 ， 本 文 使 用 了 LSMR 算法 来 解 方程 式 (14)，LSMR 算法 是 
FONG 等 人 "5 在 2010 年 提出 的 一 种 解 线性 方程 组 的 迭代 算法 。 和 LSQR 算法 一 样 ，LSMR 算法 也 是 用 于 
求 线性 方程 组 4x = b 的 解 或 者 使 Ax- 中 最 小 的 解 ， 但 和 LSQR 算法 相 比 ， 使 用 LSMR 算法 可 以 更 快 和 
更 安全 的 达到 收敛 条 件 。 
2.3 改进 的 局 部 线性 回归 估计 器 
非 参 数 SSB 模型 最 终 得 到 的 是 一 个 关于 风速 U 和 有 效 波 高 SWH 的 2 维 查 找 表 ， 需 要 计算 查找 表 中 每 
个 网 格 点 的 输入 x = (U0,SWH) 所 对 应 的 SSB 估计 值 。 虽 然 利 用 LSQR 和 LSMR 算法 来 解 线 性 方程 式 (14) 可 
以 避免 式 (15) 式 矩阵 求 道 运算 ， 但 式 (10) 中 的 加 权 函 数 avG - x,) 和 式 (19) 中 加 权 和 矩阵 A 中 的 元 素 
aw (xX, - xX») 都 是 通过 局 部 线性 回归 估计 器 得 到 , 从 式 (5) 和 式 (8) 也 可 以 看 出 通过 局 部 线性 回归 估计 器 得 到 加 
权 函 数 时 也 需要 用 到 高 维 矩阵 运算 ， 本 文 希 望 避免 高 维 矩阵 运算 ， 在 对 结果 影响 不 大 的 条 件 下 ， 能 大 大 降 
低 运 算 时 间 。 
对 于 传统 的 利用 风速 (U) 和 有 效 波 高 (SWH) 来 建立 SSB 估计 模型 时 ，x = Co,z)= (U,SWH) ，z 在 x 的 
邻 域 ， 回 归 函 数 > 的 局 部 线性 模型 可 以 表示 为 : 


r(z)= pu 下 PC 21) 十 b,(z, X,) (16) 


通过 找到 bp = (b,b,b,)1 ， 使 得 式 (17) 的 函数 最 小 : 


a 


N 2 
2=a b,- b- bo x)- bt x) JK a -x) (17) 


0 TL 


令 了 到 - 0, 到 - 0， 径 0, 那么 可 以 得 到 一 个 方程 组 : 
po Tp, $b, 
Sibot+ Sb,+ Sb, = 4 
Sibo t+ Sibit+ Sb, = 人 (18) 
Sbot+ Sib,+ Sb, = 2 


， 9， b, C 和 9 可 以 表示 为 : 
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y,Kn(X- xX,) 
je 


N 
a 
N 
b= a yx - XK ,x- x,) 
X 
a 


Ss y iXis - XK (x -Xx,) 
区 m 
S mn = d (x 一 x1) (x, 可 x,)K, (x = xX,) 
i=1 
510 So 1 曾 
520 Sn  ， 这 里 用 伪 逆 来 求 b : 


人 = pinv(S )? (20) 


际 
加 


$= ein 聂 inv(CS) 笑 


(21) 
令 Q = en ?pinv(S)， 那 么 

b=QDa QW Qe 

=a 篇 D+QOW z+ x) Kx) 

所 以 ， 权 值 函 数 可 以 表示 为 : 

a (x, -x)= 答 D+2O)?o， x)+OQO3) NU, x,) K,(x, -x) (23) 
将 式 (23) 蔡 换 到 式 (8)， 就 可 以 将 ILLR 估计 器 用 于 非 参数 SSB 模型 的 构建 ， 和 LLR 估计 器 相 比 ，ILLR 估计 
器 避免 了 高 维和 矩阵 运算 。 
昌 式 (8) 和 式 (98) 可 知 ， 两 种 估计 器 得 到 的 权 函 数 有 所 区 别 ，ILLR 估计 器 和 LLR 估计 器 在 本 质 上 是 相同 
的 ， 都 是 通过 对 式 (3) 求 偏 导 ， 并 使 其 偏 导 为 零 得 到 ， 如 果 将 式 (5) 重 写 为 下" 古 下 ,》 = 革 了 Wy ， 其 展开 后 就 
得 到 了 式 (18)， 所 以 只 要 枣 ) 印 瑟 ， 不 是 奇异 矩阵 ，ILLR 估计 器 和 LLR 估计 器 得 到 的 结果 就 完全 相同 。 在 
LLR 估计 器 中 ， 由 于 W 是 对 角 和 矩阵 ， 马 , 卫 和 半 了 W 在 运算 过 程 中 会 产生 很 多 结果 为 零 值 ,它们 会 继续 参与 
竺 WX ,和 团 5Wy 的 先 阵 运算 ， 而 在 ILLR 估计 器 中 ，LLR 估计 器 中 对 WW 竺 ,和 处 Wy 在 计算 过 程 中 产生 
的 那些 中 间 结 果 为 零 的 值 不 会 参与 式 (18)， 式 (9) 的 运算 ， 这 惑 大 大 节约 了 计算 时 间 。 从 算法 时 间 复 杂 上 度 方 面 
分 析 ， 在 LLR 估计 器 中 ， 权 函数 avG, -x) 由 T = er (XjWX,) ' 对 了?W 得 到 ，ew 是 (1”3) 的 矩阵 ， 基 ,是 


三 


(WW “3) 的 窍 阵 ,本 是 (NWN “NN) 的 和 矩阵， 处 7?WW 是 (3 N) 的 向 


相 乘 需要 3 


Wa 


N? 次 乘法 运算 ， 


(XX JW,)" 和 对 >W 相 乘 需要 3 次 乘法 运算 ， 


2N 和 (m + n+ 有 D?N 次 乘法 运 入 


OWN)。 
3 数值 实 


验 


3.1 数据 来 源 和 处 理 


本 文 使 用 Jason-2 高 度 计 的 地 球 物理 数据 集 (GDR) 产 品 ， 
TO POSDIOONS Dow MW De 被 公 
SSB 模型 时 使 用 了 cycle 51~100 
用 了 2009(cycle19~cycle55)， 

建立 SSB 模型 需要 用 到 风速 、 
需 的 各 种 校正 项 都 需要 插值 到 交叉 点 处 ， 本 文采 用 3 次 样 条 所 
在 交叉 点 两 侧 各 有 4 个 连续 的 测量 


3.2 核 函数 和 带宽 的 选择 


核 函 数 和 带宽 的 选择 对 于 非 参 数 SSB 模型 
于 非 参数 SSB 模型 的 带宽 包括 全 局 带宽 和 局 
内 速 的 带 宽 有 几 和 有 效 ; 皮 高 的 带 市 殉 he 都 为 


Epanechnikov 核 函 数 ， 而 
有 的 输入 x = (U,SWH)， 
根据 数据 密度 

和 高 斯 核 函 
了 计算 效率 ， 但 在 数 志 
从 而 在 数据 稀 疏 区 域 优 势 较 为 明 


调节 和 带 中 宽 ， 


ha 


宽 的 纪 


人 
日 


3.3 实验 


) 2 
和 局 部 带宽 ， 对 应 
结果 


证 


动 进行 调节 。 


，S 是 (3 3) 的 和 矩阵， 所 以 利用 ILLR 估 


共 50 个 cycle 的 数据 ， 而 在 评估 ILLR 估计 
ai 2011(cycle93~cycle128) 


有 效 波 高 、 交 叉 点 处 未 经 SSB 校正 的 海 


由 式 (18) 可 知 , 计算 a, bc 和 9 ，， 


E 阵 ， 轩 WX, 是 (3”3) 的 矩阵 ， 所 以 ”和 W 
和 了 和 对 , 相 肝 需 要 3?N 次 乘法 运算 ， 和 矩阵 对 WX, 求 首 需要 3 次 乘法 运算 ， 
er 和 (了 ?WW 于,) 对,W 相 乘 需要 
用 LLR 估计 器 获得 权 函 数 的 时 间 复 杂 度 为 OOV?) 。 


3N 次 乘法 运算 ， 所 以 利 
分 别 需 要 进行 N，2N， 
权 函 数 的 时 间 复 杂 度 为 


计 器 获得 


二 
年 


Jason-2 卫 


于 2008 年 6 


月 发 射 ， 它 是 


\ 认 为 目前 精度 最 高 的 测 高 J 


DD 星 。 本 文 在 建立 非 参 数 


| 器 对 SSB 估计 结果 的 影响 时 使 
£3 年 的 数据 。 


用 高 度 不 符 值 。 
上 值 的 方法 ， 为 了 得 到 更 好 的 插值 效果 ， 只 有 
点 时 得 到 的 交叉 点 才 有 效 ， 这 样 可 以 保证 较 好 的 高 度 计数 据 质量 。 


计算 海面 高 度 所 


I 非常 重要 。 目 前 ， 


最 常用 的 核 函 数 有 高 斯 核 函数 和 
部 带宽 两 种 ， 


全 局 带宽 是 指 对 于 所 


= 
Mi 


司 定 值 ， 而 局 


数 相 比 ，Epanechnikov 核 函数 可 以 将 线性 系统 式 (14) 中 和 矩 阵 B, 转换 稀 玻 
四 稀 朴 区 域 估 计 的 效果 比较 差 。 和 全 局 带宽 相 比 ， 局 部 带宽 可 以 根据 数据 的 密集 程度 
本 文 的 数值 实验 只 是 为 了 评估 改进 的 局 部 线性 估计 器 的 计 


部 带宽 是 指 h, 和 及 ， 会 


矩阵， 从 而 大 大 提高 


效率 以 及 对 估计 结果 的 影响 ， 所 以 对 核 函数 和 带宽 的 选择 没有 太 多 的 要 求 。 本 文采 


第 


1 种 是 采用 高 斯 核 函数 和 全 局 带 
于 文献 [12] 的 模型 。 


3.3.1 LLR 估计 器 和 ILLR 估计 器 估计 SSB 效率 比较 


站 -证 
A 


比较 利 


本 总 
ni 


所 用 的 时 
高 计 和 
要 先 确 
的 增 大 ， 
辣 增幅 


和 


应 的 SSB 估计 值 ， 


是 一 


的 计算 时 间 也 会 明 ! 
随 着 N 的 增 大 ff 


相对 较 小 。 
计 SSB 所 需 的 时 间 ，N 越 大 ， 改 进 的 效果 越 明 

非 参数 SSB 模型 最 终 得 到 的 是 一 个 关于 U 和 SWH 的 查找 表 ， 
通过 U 和 SWH 只 能 得 到 有 限 的 SSB 的 
ee 个 新 的 研究 方向 
高 维 查 找 表 ， 其 所 对 应 的 网 
1 并 且 随 着 维 
对 于 LLR 估计 器 是 


上 用 LLR 估计 器 估计 SSB 所 需 
从 表 1 和 表 2 可 以 看 出 ， 


[| 


巨大 的 挑战 ， 


格 点 数 将 远 
度 的 增加 ， 必 须要 


] LLR 估计 器 和 ILLR 估计 器 估计 SSB 所 月 
得 到 的 结果 ， 随 着 每 次 用 于 SSB 估计 的 交叉 点 数 N 的 增 大 ，LLR 估计 器 有 
显 增 大 。ILLR 估计 器 中 和 矩阵 的 维 
而 增 大 , 但 其 增 大 的 幅度 相对 较 小 。 表 2 
间 ， 正 常情 况 下 ，Epanechnikov 核 函 数 可 以 将 线性 系统 式 (14) 中 的 入 


了 两 种 核 函 数 和 带 


带宽 ， 对 应 于 文献 [10] 的 模型 ， 第 2 种 是 Epanechnikov 核 函数 


的 时 间 ， 表 1 给 出 了 利用 高 斯 函数 和 全 局 带 
的 矩阵 的 维 


度 就 相应 地 增 大 ， 所 


度 不 会 随 着 N 的 增 大 而 改变 ， 


效率 ， 降 低估 计 SSB 所 需 的 时 间 ， 但 由 于 采 
定 其 所 在 网 格 所 对 应 的 带宽 ， 这 样 就 会 增加 计算 量 ， 从 而 使 得 计算 时 
的 计算 时 间 
和 LLR 估计 器 相 比 ， 利 ) 


显 。 


[6, 16, 17] 。 


用 了 占 计 SSB 时 ， 


间 增 
ILLR 估 


局 部 人 带宽， 在 


大 幅 增 加 ， 而 利用 


虽然 其 所 需 的 时 间 也 会 


给 出 了 利用 Epanechnikov 核 函 数 和 局 部 带宽 估计 SSB 
E 阵 B, 转换 为 稀 疏 和 矩阵， 从 而 提 


对 于 每 一 输入 x ， 都 需 


大 。 和 表 1 一 样 ， 随 着 N 
计 器 估计 SSB 所 需 的 计算 时 
jj ILLR 估计 器 大 可 以 大 幅度 地 降低 估 


言 息 ， 将 新 的 参量 (如 波 


4 的 构建 才 色 


这 一 问题 


表 1 利用 高 斯 核 函 数 和 全 局 带宽 估计 SSB 所 用 的 时 间 


估计 器 
LLR 
ILLR 


174 S 
28 S 


N=500 


N=2000 
37 min 
1.5 min 


需要 计算 每 个 查 


找 表 中 每 个 网 格 点 所 对 
周期 参数 ) 引 入 SSB 模型 


在 建立 高 维 非 参数 SSB 模型 时 ， 最 终 得 到 的 结果 将 
远 多 于 2 维 查 找 表 的 网 格 点 数 ， 使 用 LLR 估计 器 所 需 消耗 的 
有 更 多 的 数据 参与 模型 
而 ILLR 估计 器 则 可 以 很 好 地 解决 ; 


保证 模型 的 有 效 性 ， 这 


N=5000 


6.8h 
6 min 


表 2 利用 Epanechnikov 核 函 数 和 局 部 带宽 估计 SSB 所 用 的 时 间 


估计 器 
LLR 
ILLR 


N=500 
330 S 
38S 


3.3.2 LLR 估计 器 和 ILLR 估计 器 估计 SSB 结果 比较 


图 1 是 利用 高 斯 核 函数 和 全 局 
部 分 表示 高 度 计 测 
1(a) 和 图 1(b) 的 等 值 线 分 
等 值 线 的 单位 为 cm。 很 明显 , 利用 


大 约 在 10-+cm 
图 
估计 器 和 ILLR 
两 种 估计 器 得 到 


数 得 到 的 结果 平滑 。 图 2 


到 的 SSB 基本 相同 。 
图 3 给 出 了 利用 LLR 信 t 


佛 宽 得 到 的 结果 ， 


别 表示 利用 LLR 估计 器 和 
两 种 估计 器 得 到 
条 件 下 ，SSB 的 幅 值 都 随 着 SWH 的 增 大 而 增 大 ; 
是 先 增 大 ， 而 当 U 增 大 到 一 定 值 时 ，SSB 的 幅 值 都 逐渐 减 小 。 
汕 计 器 得 到 的 SSB 之 差 在 (U,SWH) 平 面 的 分 布 ， 可 以 看 出 ， 
量 级 ， 基 本 上 可 以 忽略 。 

2 是 利用 Epanechnikov 核 函数 和 局 部 带宽 得 


十 器 和 ILLR 估计 器 得 到 的 SSB 估 


计 器 得 到 的 SSB 查找 表 分 别 用 于 2011 年 Jason-2 高 度 计 


估 1 致 性 
估计 ， 


非常 好 。 
分 别 利用 每 年 的 数据 进行 


和 局 部 带宽 
不 管 是 哪 种 核 函 数 和 带宽 
目前 ， 


的 组 合 ， 


统计 ， 表 3 给 
器 得 到 的 SSB 估计 值 之 差 的 最 大 值 在 10-3 cm 是 


再 将 两 种 估计 器 得 到 的 的 查找 表 分 别 月 


值 的 方差 减 去 经 过 SSB 校正 
能 够 解释 的 部 分 ， 解 释 方 差 越 大 ， 

2009~2011 3 年 的 数据 计算 1 用 两 种 
的 解释 方差 之 差 。 当 
而 当 使 用 Epanechnikov 核 函数 和 局 部 
所 以 不 管 是 哪 种 核 函数 和 带宽 的 组 合 


当 使 用 高 斯 核 函 数 和 全 


根据 数据 的 密集 程度 将 数 
量 值 个 数 小 于 100 的 区 域 ， 代 表 数 据 稀疏 区 域 ， 而 另 外 一 部 分 则 代表 数据 密集 区 域 。 图 
ILLR 估计 器 得 到 的 SSB 查找 表 在 
的 SSB 随 U 和 SWH 的 变化 趋势 完全 相同 , 在 U 相同 的 
而 在 SWH 相同 的 条 件 下 ， 随 着 U 


到 的 结果 ， 图 2(a) 和 
估计 器 得 到 的 SSB 查找 表 在 (USWH) 平 
的 SSB 随 U 和 SWH 的 变化 趋势 也 是 相 
高 海 况 和 高 风速 低 海 况 的 数据 稀疏 区 域 ， 由 于 使 用 到 的 数据 比较 少 ， 
2(c) 的 等 值 线 表示 利 月 


N=2000 N=5000 
12.4h 


8.9 min 


77 min 
2 min 


的 
图 1(c) 的 等 值 线 表示 利用 
在 数据 密集 区 域 ， 两 种 估计 


| 器 得 


的 分 布 ， 等 值 线 的 单位 为 cm， 同 图 


四 划分 为 两 部 分 


U,SWH) 平 


得 到 的 结果 不 如 


图 1 中 利 ) 
日 LLR 估计 器 和 ILLR 估计 器 得 到 的 SSB 之 差 在 (U,SWH) 平 
分 布 ， 在 数据 密集 区 域 ， 利 用 两 种 估计 器 得 到 的 SSB 之 差 大 约 在 10-6cm 量 级 ， 可 以 认为 利用 两 也 


， 图 中 阴影 


用 的 分 布 ， 


曾 大 ，SSB 的 幅 值 都 
LLR 估计 器 和 ILLR 
到 的 SSB 之 差 


图 2(b) 的 等 值 线 分 别 表示 利用 LLR 
1 一 样 ， 利 用 
同 的 ， 由 于 使 用 了 Epanechnikov 核 函 数 ， 在 低 


风速 
高 斯 核 函 
的 
计 器 得 


估 


计 值 的 散 点 图 ， 


图 中 的 结果 是 通过 将 两 种 估 


SSB 的 估计 得 到 , 可 以 看 出 , 两 生 


的 影响 都 可 以 忽略 。 


估计 器 得 
有 于 2009~2011 年 Jason-2 高 度 计 SSB 的 
出 了 统计 结果 ， 当 利用 高 斯 核 函 数 和 全 局 带宽 
量 级 ， 平均 值 在 10-* cm 寺 
时 ， 两 种 估计 器 得 到 的 SSB 估计 值 之 差 的 最 大 值 在 10-5 
ILLR 估计 器 对 SSB 的 估计 结果 


时 ， 两 种 估 
量 级 ， 而 当 利 用 Epanechnikov 核 函 数 
cm 量 级 ， 平 均值 在 10-6 cm 量 级 。 所 以 


到 的 SSB 


证 


里 解 为 海 
同 的 核 函 数 和 带 


赴 计 器 得 到 的 SSB 


忽略 ， 从 而 也 不 会 影响 高 度 计 最 终 的 测 高 精度 。 


带宽 时 ,两 种 估计 器 得 至 
， 和 LLR 估计 器 相 比 ， 


的 解释 方差 , 表 4 给 出 了 利 | 


ILLR 估计 器 对 SSB 模型 


到 际 上 主要 采用 解释 方差 来 对 SSB 模型 进行 评估 。 解 释 方差 是 指 未 经 SSB 校正 的 海面 高 度 不 符 
后 的 海面 高 度 不 符 值 的 方差 , 解释 方差 可 以 
说 明 SSB 模型 越 有 效 。 对 两 种 不 


高度 不 符 值 的 方差 中 SSB 
宽 组 合 ， 
用 两 种 估计 器 得 到 的 SSB 
局 带宽 时 , 两 种 估计 器 得 到 的 SSB 的 解释 方差 之 差 在 10scms 量 级 ， 

| 的 SSB 的 解释 方差 之 差 在 10-9~10- cm? 
J 有效 性 的 影响 都 可 以 


分 别 使 用 


量 级 ， 
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(c) 两 种 估计 器 估计 的 SSB 之 差 (105cm) 


日 Epanechnikov 核 函 数 和 局 部 带宽 得 到 的 结果 在 (U,SWH) 平 面 的 分 布 
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(b) 利 用 Epanechnikov 核 函数 和 局 部 带宽 得 到 的 结果 


日 LLR 估计 器 和 ILLR 估计 器 得 到 的 SSB 估计 值 的 散 点 图 
表 3 利用 两 种 估计 器 得 到 的 SSB 之 差 (cm) 统 计 


2009 年 2010 年 2011 年 2009 年 2010 年 2011 年 
高 斯 核 + 全 局 带宽 1.26X10- 1.01X10-? 1.13X10- 5.45X10- 5.51X10- 5.49X10- 
Epanechnikov 核 + 局 部 带宽 3.63X1075 8.19X1075 5.01X1075 2.39X10-6 2.41X10-6 2.40X10-6 
P 核 可 外 市 殉 
表 4 利用 两 种 估计 器 得 到 的 SSB 的 解释 方差 之 差 (cm?) 
检验 数据 集 
核 函 数 和 带宽 组 合 一 — 
2009 年 2010 年 2011 年 
高 斯 核 + 全 局 带宽 4 和.83X10-5 4.88X10-5 5.27X10-5 
Epanechnikov 核 + 局 部 带宽 7.88X10-8 1.28X107 1.77X10”9 


结论 


时 ， 


非 参 数 模型 由 于 其 较 高 的 精度 已 被 广泛 用 于 雷达 高 度 计 的 SSB 校正 ， 在 利用 非 参 数 估计 方法 估计 SSB 
通常 会 采用 LLR 估计 器 来 求 加 权 和 矩阵 里 的 权 函 数 ， 而 LLR 估计 器 涉及 高 维和 矩阵 运算 ， 当 数据 量 较 大 时 


需要 大 量 的 计算 时 间 。 


使 月 
利 月 


本 文 提出 了 一 种 ILLR 估计 器 ， 该 估计 器 可 以 避 免 高 维 算 阵 运算 ， 从 而 极 大 地 提高 计算 效率 。 本 文 分 别 
高 斯 核 函数 加 全 局 带宽 以 及 Epanechnikov 核 函数 加 局 部 带宽 这 两 种 组 合 来 估计 SSB， 在 两 种 组 合 中 ， 
日 ILLR 估计 器 得 到 的 SSB 和 利用 LLR 估计 器 得 到 的 SSB 之 差 基 本 上 可 以 包 略 ;，ILLR 估计 器 可 以 大 幅 


度 地 降低 估计 SSB 所 需要 的 时 间 ， 当 数据 量 很 大 时 ， 这 种 改进 的 效果 更 为 明显 。 目 前 业务 化 运行 的 高 度 计 
都 是 采用 基于 风速 和 有 效 波 高 的 2 维 SSB 模型 ， 将 新 的 参量 引入 SSB 模型 已 经 成 为 SSB 经 验 模型 研究 的 一 
个 新 方向 ， 本 文 提出 的 ILLR 估计 器 可 以 在 不 影响 SSB 估计 结果 条 件 下 大 幅度 降低 估计 SSB 所 需要 的 时 间 ， 


这 为 实现 高 维 非 参 数 SSB 模型 的 实时 运算 葛 定 了 基础 。 
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究 方 向 为 雷达 高 度 计数 据 处 理 . 
士 生 导 师 ， 主 要 研究 方向 为 星 载 雷达 高 度 计 系统 技术 、 合 成 孔径 雷达 高 度 计 系统 技术 和 信号 处 理 
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